Role of relaxation in spin Hall effect 
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The role of the relaxation due to the impurity scattering and/or the contacts to leads/electrodes 
are studied for the spin Hall effect (SHE). Relaxation is essential to attain the steady state and 
also to the spin accumulation, but has been considered to be harmful for the intrinsic SHE (1SHE). 
These issues are examined quantitatively on two types of 2D models, i.e., (a) Rashba model for 
n-type GaAs, and (b) Luttinger model for p-type GaAs. It is found that ISHE is robust against 
the realistic strength of the disorder producing the observable amount of the spin accumulation. 
Especially in model (b) the spin current and the accumulation are an order of magnitude larger 
than those in model (a). Experimental observations are discussed quantitatively from these results. 
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The spin Hall effect (SHE) is a new realm of spintron- 
ics, by which spin current is produced perpendicular to 
the applied electric field. This enables the spin injection 
to the semiconductors without the magnets or magnetic 
field. This possibility has been proposed long ago QUEL 
but recent intensive interest comes from the theoretical 
proposal that the intrinsic mechanism due to the topo- 
logical nature of the Bloch wavefunction in the presence 
of the spin-orbit interaction (SOI) gives rise to orders 
of magnitude larger effect in conventional semiconduc- 
tors such as GaAs 0, 0, 0]- Recent two papers 0, @| 
reported on the experimental observation of SHE. One 
is on n-type GaAs where the charge current J produces 
the spin accumulation detected by Kerr rotation spec- 
troscopy near the edges of the sample transverse to J • 
These authors concluded that the effect is due to the 
extrinsic origin since the effect is rather independent of 
the orientation of the sample. The other is on p-type 
GaAs where also the spin accumulation is detected by 
the circularly polarized LED || . Their sample has more 
carrier density compared with above, and the estimation 
of the relaxation is small which lead the authors to con- 
clude the intrinsic origin. However the debates on the ori- 
gin of the SHE, i.e., the intrinsic or extrinsic, continues. 
This situation is in parallel to the anomalous Hall effect 
(AHE), where the long standing controversy between the 
extrinsic impurity induced mechanism (such as the skew 
and side jump model [11|) and intrinsic 
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one [l2l Il3l Il4l Il5| still continues. Therefore it is of vi- 
tal importance to study the effect of impurity scatterings 
and/or relaxation on the SHE quantitatively taking into 
account the realistic values of parameters and experimen- 
tal setups, which we undertake in this paper. Actually 
there are several prece ding wo rks address ing the issue of 
relaxation in SHE [H EjEl E |H HJ H H |H . 
However many of the works focus on the limiting case 
of weak disorder El [H H [3 |H H EL or non- 
equilibrium state is not taken into account |23l Eif , or 



lacking a quantitative comparison with existing experi- 
ments |25|. Also the studies on the 4-band model which 
describes p-type GaAs are missing except ref. 20] . This 
is partly due to the fact that the spin relaxation is so 
rapid for p-holes since the SOI is much larger at the top 
of the valence bands compared with the bottom of the 
conduction bands, which has been assumed to be dis- 
advantageous for the spin accumulation. However, the 
strong SOI promotes the spin current and in any case 
the relaxation is needed to produce the spin accumula- 
tion, which is time reversal (T)-odd, from the T-even spin 
current. Therefore it is a nontrivial issue which is more 
advantageous, p-type or n-type GaAs. 

To answer all these questions, we study in this paper 
the Rashba and Luttinger models defined on the square 
lattice in terms of the Keldysh formalism applied to the 
finite size sample attached to the leads/electrodes jH 
[27]]. These models reproduce the continuum version of 
each system near the T-point of Brillouin zone, which 
is most relevant for the low carrier density. The energy 
unit t (hopping parameter) and the length unit a (lattice 
constant) will be fixed later when we compare our re- 
sults with experimental data. The model for the Rashba 
system is expressed as 
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When the Fermi energy is near the band bottom and 
S <C 1, the effective mass m* and the Rashba coupling a 
are given by m* ~ l/(2ta 2 ) and a ~ Sta, respectively. 
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The model for the Luttinger system is defined by 
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where To is the 4x4 unit matrix and other T-matrices 
are defined in ref. 0. Assuming p-type GaAs thin layer, 
we shall take S — 0.29, M = 2t. This parameter 
set is corresponding to the typical Luttinger parameters 
7l : 72 : 73 = 6.92 : 2.1 : 2.1 for GaAs and ((k z a) 2 ) = 
M I (2St) ~ 1.86, is determined by the profile of the con- 
fined wavefunction along the z-direction |29l | . In the origi- 
nal three-dimensional system, the kinetic terms contains 
Ti^-matrices and t r ' r , ^ 0. We can approximately ne- 
glect these terms in the quasi-two-dimensional system 
confined in a thin layer |29|. When i ' , = 0, the 4x4 
T-matrix space will be decoupled to two of the 2x2 ma- 
trix space. However, as long as we consider the doped 
system, this point does not lead to any crucial difference. 
We shall take the units in which H = c= 1. 

For each model, we obtain the Keldysh matrix Green 
function by solving the integral equations numeri cally 
in the self-consistent Born approximation |2rl \2l\ . 
The retarded self-energy is given by SLir^^^) = 

5 ^j«MK'1 1 where Z R > cont (E) is the 
contact self-energy, and r r (E) the local lifetime due 
to disorders. The local lifetime r r (E) is determined 
self-consistently by the recursion equation, - ^ = 
7 N r (E) , where 7 represents the strength of disorder and 
N r (E) = jfaTrW[G*.(E) ~ G$ r (E)] is the local den- 
sity of states per unit cell. The lesser Green function 
G< r (E) is also determined self-consistently, by which the 
spatial-dependent physical quantities, i.e, the spin den- 
sity/current, and charge density/current, can be calcu- 
lated. We take the sample of finite size L x x L y with 
the electrodes attached at x = ±L x /2, while the open 
boundary condition is imposed in the y-direction [27| • It 
is noted here that the edge modes do not play crucial role 
in the doped case. We take the small chemical potential 
difference dfi/L x = 5 x lQ~H/a to study the linear re- 
sponse regime. The chemical potential in equilibrium is 
taken as fj,o = — 3.5t for the Rashba system and [1q = for 
the Luttinger system. Then, the Fermi energy mea- 
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sured from the band edge is ~ 0.5i for the Rashba 
system and ~ 2t for the Luttinger system respec- 
tively. In both systems, the dispersion near the Fermi 
level is almost quadratic and isotropic. In the Luttiger 
system, the Fermi level is crossing only the heavy hole 
bands. It is noted that our models describes only the 
electronic states near the T-point in Brillouin zone but 
not those of whole Brillouin zone. Therefore, the above 
values of do not necessarily mean large carrier con- 
centrations. The substantial amount of carrier concen- 
tration is determined after fixing the parameters t and 
a, which will be done when we compare our results with 
experimental data. 
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FIG. 1: Spin accumulation 5S Z and the spin current SJ S " in 
(a), (b) the Rashba system (S=0.05) with 7 = 2 (1/r ~ 0.3t), 
and the Luttinger system with (c), (d) 7 = 1 (1/r ~ 0.25t) 
and (e), (f) 7 = 5 (1/r ~ 1.25t). The system size is L x x 
L y — 40a x 40a for the Rashba system and 40a x 20a for 
the Luttinger system. The chemical potential of electrons 
is jito + S/J./2 at x = —L x /2 and /io — Sfi/2 at x = L x /2. 
The charge current 8 J (not shown) flows in the negative x- 
direction. 

We show the obtained results in Figs. 1,2, and 3. Fig.Q] 
shows the accumulation pattern of spin z-component 
SS Z and the spin current SJ S * for (a),(b) the Rashba 
(S = 0.05) and (c),(d),(e),(f) Luttinger systems. Here 
we take the definition of the spin current as J r £, = 
\ (S^Jrr' + Jrr'S/i), where S is the spin-^ matrices for 
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the Rashba system and the spin-| matrices for the Lut- 
tigner system, and J rr ' is the charge current. The disor- 
der strength is taken as 7 = 2 for the Rashba system, and 
7=1 and 5 for the Luttinger system. Then, the inverse 
lifetime is 1/r ~ 0.3t for the Rashba system, 1/r ~ 0.25i 
and 1.25i for the Luttinger system. Figure [21 shows (a) 
the spin z-component 5S Zl (b) the charge current 5J X 
and (c) the divergence of the spin current V • S J z at the 
x = cross-section, while Fig. 13 (a) the spin y-component 
SS y , (b) the spin current SJy and (c) the electron den- 
sity 5n at the y = cross-section. 
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FIG. 2: Distribution of (a) 5S Z , (b) 5J X and (c) V • 5J Sz at 
x — in the Rashba and Luttinger systems. The system size 
is L x x L y = 40a x 40a for the Rashba system and 20a x 40a 
for the Luttinger system. Only the region y < is shown. 
SS Z and V ■ SJ' Sz are the odd functions of y. The disorder 
strength is 7 = 2 for the Rashba system and 7=1 (solid) and 
7 = 5 (dashed) for the Luttinger system, respectively. As for 
the Luttinger system, 5S Z in 7 = 1 and V -8J S * in 7 = 1 and 
5 are rescaled. 

As shown in Fig.^( a ), the spin accumulation occurs all 
along the edges for Rashba system. This is in sharp con- 
trast to the theoretical prediction for the Rashba system 
in ref . |2l| , where the spin current and spin accumulation 
is finite only near the electrodes. However since the spin 
current is not the conserved quantity, not all the spin 
current contributes to the spin accumulation. There ap- 
pears sink (V • SJ S * < 0), which is related to the spin 
torque density and the spin relaxation [3(| j near the neg- 
ative y-edge , where the spin current is partly absorbed 
as shown in Fig. [21 (c). The spin accumulation comes 
from the remaining part of V • 5J S - which is not can- 
celed by the torque density and balancing with the spin 
relaxation. Therefore, the bulk SHE, the direction of the 
spin current and the sign of the spin accumulation are 
consistent. 

As for the Luttinger system, the spin current is 
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FIG. 3: Distribution of (a) 6S y , (b) 5J%* and (c) Sn at y = 
in the Rashba and Luttinger systems. The system size and 
the disorder strengths are the same as in Fig. H SJS* of the 
Luttinger system is rescaled both in the cases of 7 = 1 and 5. 
Only the region x < is shown. Sn is the odd functions of x. 



strongly enhanced near the electrodes x = ±L x /2, but 
the resultant spin accumulation is rather flat along the 
edges y = ±L y /2. The sign of spin accumulation pattern 
for 1/r < t (7 = 1) is opposite to what is expected from 
that of spin Hall conductivity in the bulk while the di- 
rection of spin current in the bulk is consistent. This is 
because the spin current near the contact gives the op- 
posite contribution in this system as shown by the pro- 
file of 8Jy z in Fig. [21 (b). When the relaxation effect is 
increased, the bulk property becomes dominant even in 
a small system and the accumulation pattern coincides 
with what is expected from the spin Hall conductivity. 
We need more rigorous argument on the definition of 
conserved spin current |3l| in order to investigate this 
problem furthermore. As for the Rashba system, as seen 
in Fig. [3| (a), the in-plane spin accumulation 6S y perpen- 
dicular to the electric field is finite in the bulk as dis- 
cussed in refs. [32I l33| . On the other hand, there appears 
no in-plane spin accumulation for the Luttigner system 
as long as t r ' r , = 0, i.e., no hybridization between the 
decoupled 2x2 matrix spaces in 4 x 4 T-matrix space. 

Next, Fig. 01 shows the inverse lifetime dependence of 
5S Z . The disorder strength is taken as 7 = 0.05, 0.1, 0.2, 
0.5, 1, 2, 5 and 10. Sample points are the peak values of 
\SS Z \ in the region, |x| < 2a and \y\ ~ L y /2. It is seen 
that the spin accumulation is larger for stronger SOI. 
This means that the magnitude of the spin current is the 
more important factor than the spin lifetime. Therefore, 
it is concluded that p-type GaAs is more advantageous 
than n-type to observe the spin accumulation due to the 
ISHE. Another observation is that the spin accumulation 
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is rather robust against the relaxation up to 1/r ~ O.lt. 

S|i/L v = 5.0x Iff 4 tla 
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FIG. 4: Spin accumulation as a function of the inverse life- 
time. Sample points are the peak values of \SS Z \ in the region, 



lad < 2a and 



Ly/2. The colors are corresponding with 



those in Figs. |5]and [3] The sign of the accumulation pattern 
in the Luttinger system changes around the dashed vertical 
line. 

In order to compare the results with the experiment |2J 
on n-type semiconductor, we fix t and a in the Rashba 
system as t — 4 meV and a = 14 nm. The effective 
mass and the carrier density is estimated as m* ~ 0.05m e 
and «3D ~ 3 x 10 16 cm -3 , respectively. The Rashba 
spin splitting is Ar ~ 5.75 meV (S <C 1). The applied 
electric field in our simulation is E ~ 0.14 mV/im _1 . In 
the case with 7 = 2 which corresponds to 1/r ~ 0.3i, The 
charge resistivity is p c ~ 144 f2/im, and the spin Hall 
resistivity multiplied by 1/e is ~ 2.9 x 10 5 , 6.4 x 
10 4 , 1.2 x 10 4 and 4.1 x 10 3 fi/im for S = 0.01, 0.02, 
0.05 and 0.10 respectively. When the results are linearly 
extrapolated to E = 10 mVjum" 1 , we obtain : \SJ X \ ~ 70 
pApm' 2 , \e6J*f*\ ~ 35, 160, 870, 2400 nA^m" 2 , and 
\SS Z \ ~ 2.6, 10, 64, 150 /1m- 3 for S = 0.01, 0.02, 0.05 
and 0.1, respectively. Although the Rashba coupling for 
S = 0.01 is still an order of magnitude larger than the 
experimental estimation of the Rashba coupling in ref.0, 
Pc 5J X and SS Z are consistently corresponding to the 
experimental data. However, p s and 8Jy z are not. This 
may be attributed to the spin source and sink where V • 
SJ Sz ^ 0. In the Luttinger system, t and a are fixed as 
t = 10 meV and a = 5 nm, which correspond to the case 
with ~ 20 meV and n 2 D ~ 1.6 x 10 12 cm~ 2 . In the 
case with 7=1, the charge resistivity is p c ~ 18 U/im. 
In the experiment a current I p ~ 100 /iA is applied 
to the p-channel of 1.5 /im width. Although we cannot 
exactly read the depth of the channel from the reference, 
it is considered to be of the order of the length unit in the 
model, i.e., a = 5 nm. In our system, this corresponds 
to E x ~ 200 mV/zru -1 . When the results are linearly 
extrapolated to E x = 200 mV/mi -1 , |55 Z |2D ~ 300 /im -2 
and the spin polarization of the holes is about 2 % which 
is of the order of the observed circular polarization (~ 1 



%) of light emitted from the LED. 

Finally, it is worthwhile to estimate the inverse spin 
Hall effect in which a gradient of external magnetic 
field, i.e., a spin force, induces a transverse charge cur- 
rent and/or electric field |3l| . For the linearly modu- 
lated magnetic field, V y B z — const., applied on an open 
system, the induced electric field is given by l-B^I = 
\<7°yP c gfiBVyB z \, where g is the ^-factor and a^y 1S tne 
conductivity of inverse SHE which has the relation be- 
tween the spin Hall conductivity a s ^ x as a c x s y — — <t*£, via 
the Onsager's relation. In the case in which V y B z = 1 
Tern- 1 , p c = lOOf^m and l/\ea c / y \ = lOOOf^m, the in- 
duced electric field is of the order of lO/iVcm -1 , assuming 
g ~ 2. Although the above estimation is not large, we 
can enhance this inverse effect by lightly doped spin Hall 
insulators, i.e., materials with low carrier density and 
large SOI. This is because larger electric resistivity and 
spin Hall conductivity are advantageous for this effect. 

In conclusion, we have numerically investigated the 
spin accumulation due to the intrinsic spin Hall effect 
in the Rashba and Luttinger systems by using Keldysh 
formalism. The distribution of the accumulated spin and 
the charge and spin currents are obtained, which are com- 
pared with those of the recent experiments obtaining the 
quantitative agreement. 
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